The quantum critical behavior of antiferromagnetic itinerant systems with van Hove 

singularities of electronic spectrum 
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The interplay of magnetic and superconducting fluctuations in two dimensional systems with van 
Hove singularities in the electronic spectrum is considered within the functional renormalization 
group (fRG) approach. While the fRG flow has to be stoped at a certain minimal temperature 
2rg"i w6 study temperature dependence of magnetic and superconducting susceptibilities both, 
above and below T^q", which allows to obtain the resulting ground state phase diagram. Close 
to half filling the fRG approach yields two quantum phase transitions: from commensurate an- 
tiferromagnetic to incommensurate phase and from the incommensurate to paramagnetic phase, 
the region of the commensurate magnetic phase is possibly phase separated away from half filling. 
Similarly to results of Hertz-Moriya-Millis approach, the temperature dependence of the inverse 
(incommensurate) magnetic susceptibility at the quantum phase transition from incommensurate 
to paramagnetic phase is found almost linear in temperature. 
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I. INTRODUCTION 

The quantum critical points (QCP) in itinerant mag- 
nets have being investigated during long time. Moriya 
theorjii was first attempt to describe thermodynamic 
properties near QCP. This theory was further developed 
within Hertz-Millis renormalization group approach^. In 
more than two dimensions Hertz-Moriya-Millis approach 
predicts that the magnetic transition temperature Tc de- 
pends on the distance S to the QCP as ~ gz/{d.+z~2)^ 
In two-dimensional systems, the long-range magnetic or- 
der at finite temperatures is prohibited according to the 
Mermin- Wagner theorem, but the quantum phase tran- 
sition is accompanied by vanishing of the temperature of 
the crossover to the regime with exponentially large cor- 
relation length (renormalized-classical regime), T* 6. 

The applicability of Hertz-Moriya-Millis (HMM) ap- 
proach to magnetic systems was however recently 
questioned because of expected strong momenta- and 
frequency dependence of the paramagnon interaction 
vertices^ and possible non-analytical dependence of the 
magnetic susceptibility^ which arises due to strong 
electron-paramagnon interaction. Studying the problem 
of quantum critical behavior of these systems in terms 
of fermionic degrees of freedom may be helpful to obtain 
concentration dependence of the crossover temperature 
and provide valuable information about their magnetic 
properties near quantum critical points. 

Itinerant systems with van Hove singularities in the 
electronic spectrum have strong momentum dependence 
of interaction vertices due to peculiarity of the elec- 
tronic dispersion, and, therefore represent an interest- 
ing example for studying the quantum critical behavior. 
The competition of different kinds of fluctuations, and 
even long range orders is important in the presence of 
van Hove singularities, which makes formulation of ef- 
fective boson-fermion theories rather complicated. At 
the same time, fermionic approaches can treat naturally 
both, (anti)ferromagnetic and superconducting fluctua- 



tions, which were considered to be important near mag- 
netic quantum phase transitions in systems with van 
Hove singularities in the electronic spectrum. 

The simplest mean-field analysis of the Hubbard model 
is insufficient to study quantum critical behavior; due 
to locality of the Coulomb repulsion in this model it 
is also unable to investigate the range of existence of 
unconventional (e.g., d- or p-wave) superconducting or- 
der, and introduction of the nearest-neighbor interaction 
is required in this approach^. To study the competi- 
tion of magnetism and superconductivity in the Hub- 
bard model, more sophisticated approaches, e.g. cluster 
methods^ii and functional renormalization group (fRG) 
approaches2i2ii2iiii were used. The fRG approaches are 
not limited by the system (cluster) size and offer a possi- 
bility to study both, magnetic and superconducting fiuc- 
tuations, as well as their interplay at weak and interme- 
diate coupling. 

The fRG approaches were initially applied to the 
paramagnetic non-superconducting (symmetric) phase to 
study the dominant type of fluctuations in different re- 
gions of the phase diagram^'^ i^^'^"'^ . Although these ap- 
proaches suffered from the divergence of vertices and sus- 
ceptibilities at low enough temperatures near the mag- 
netic or superconducting instabilities, comparing suscep- 
tibilities with respect to different types of order at the 
lowest accessible temperature provided a possibility to 
deduce instabilities in different regions of the phase di- 
agram. In fact, the temperature where the vertices and 
susceptibilities diverge in the one-loop approach, is re- 
lated to above discussed temperature T* of the crossover 
to the 'strong-coupling' regime with exponentially large 
magnetic correlation length. 

To access the region T < T*, the combination of the 
fRG and mean-field approach was proposed in RefJ^, 
which was also able to study possible coexistence of mag- 
netic and superconducting order at T = (the mag- 
netic order parameter was assumed to be commensurate). 
More sophisticated fRG approach in the symmetry- 
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broken phasei^ was developed recently to avoid appli- 
cation of the mean-field approach after the RG flow; the 
application of this method was however so far restricted 
by the attractive Hubbard model, because of complicated 
structure of the resulting renormalization group equa- 
tions. 

So far only susceptibilities corresponding to spin and 
charge fluctuations with commensurate wavevectors, as 
well, as to superconducting fluctuations were carefully 
investigated. In the present paper we use the fRG ap- 
proach in the symmetric phaseiSiii and perform an ac- 
curate analysis of temperature dependence of suscepti- 
bilities with respect to both, commensurate and incom- 
mensurate magnetic order, as well as superconducting 
order. We propose extrapolation method which allows 
us to study thermodynamic properties both above and 
below the temperature at which the fRG flow is stopped, 
and extract the crossover temperature T*. This gives 
us a possibility to obtain phase diagram, capturing sub- 
stantial part of the fluctuations of magnetic and super- 
conducting order parameters without introducing sym- 
metry breaking. Contrary to the functional renormal- 
ization group analysis in the symmetry broken phase^^, 
the presented method can be easily generalized to study 
instabilities with different type of the order parameters. 



II. METHOD 



We consider the 2D t-t' Hubbard model Hn 
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~'U')N with 
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(1) 



where — t for nearest neighbor (nn) sites i, j, and 
tij = ~t' for next-nn sites {t, t' > 0) on a square lattice; 
for convenience we have shifted the chemical potential /i 
by it'. We employ the fRG approach for one-particle ir- 
reducible generating functional and choose temperature 
as a natural cutoff parameter as proposed in Refji?. This 
choice of cutoff allows us to account for excitations with 
momenta far from and close to the Fermi surface. Ne- 
glecting the frequency dependence of interaction vertices, 
the RG differential equation for the interaction vertex 
Vt = F(ki, k2, ka, k4) has the for] 

dVr dLpp 
dT ^ dT ^ 
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Vt o ,^ o Vt ■ 
dT 



(2) 



where o is a short notation for summations over inter- 
mediate momenta and spin, momenta are supposed to 
fulfill the momentum conservation law ki +k2 = k.3 +k4. 



lower signs for the particle-particle (ipp) bubbles, re- 
spectively. Eq. is solved with the initial condition 
(ki, k2, ks, k4) — U; the initial temperature is cho- 
sen as large as Tq = lO^i. The evolution of the vertices 
with decreasing temperature determines the temperature 
dependence of the susceptibilities according to^ 



dT 

dTZ'^ 
dT 



E-U'n-nrn dLph,pp(k', zfck' -|- Qm) 
'<-k' '<-±k'+q™ ^ > 14J 



Here the three-point vertices TZ^ describe the propa- 
gation of an electron in a static external field, m de- 
notes one of the instabilities: antiferromagnetic (AF) 
with q„i = (7r,7r), incommensurate magnetic (Q) with 
the wave vector q™ = Q, or d-wave superconducting 
(dSC) with — (upper signs and ph correspond to 
the magnetic instabilities, lower signs and pp to the su- 
perconducting instability); 



1^'^ - / "^^riK k', k' + q„,) m = AF or Q, 
' - \ VT(k, -k + q„,, k') m = dSC. 



(5) 



The initial conditions at To for Eqs. ([4]) are TZ^ — /k 
and Xm — 0, where the function /k belongs to one of 
the irreducible representations of the point group of the 
square lattice, e.g. /k = 1 for the magnetic instabili- 
ties and /k = {coskx — cos ky)/ A for the d-wave super- 
conducting instability, with a normalization coefficient 
A = (l/iV)Ek/k- To solve the Eqs. Q and g]) we 
discretize the momentum space in Np — 48 patches using 
the same patching scheme as in Ref J^. This reduces the 
integro-differential equations ([2]) and ([4]) to a set of 5824 
differential equations, which were solved numerically In 
the present paper we perform the renormalization group 
analysis down to the temperature Trq", at which vertices 
reach some maximal value (we choose Vi„ax = 18t). 

To obtain the behavior of the susceptibilities at T < 
^Rcf ^® extrapolate obtained temperature dependence 
of the inverse susceptibilities by fitting this dependence 
above (but close to) T^q" by polynomials of 5-th to 7- 
th order. We identify the crossover temperature T^ to 
the regime of strong correlations of the order parameter 
denoted by m from the condition that the extrapolated 
Xni'i'^m) ~ (we assume that the susceptibilities are al- 
most analytic functions of temperature in the crossover 
regime). We have checked that the obtained essen- 
tially depends on neither the order of polynomial, used 
for the fitting, nor on the fitting range. Studying the be- 
havior of T^ as a function of electron density, interaction 
strength etc. allows us to obtain the phase diagram. 



■^^ph,pp(k,k') = 



/T(ek)-/T(±£k') 



f^k T £k' 



(3) 



III. RESULTS 



and /t(£^) is the Fermi function. The upper signs 
in Eq. ^ stand for the particle-hole (Lph) and the 



We consider first small interaction strength U = 2.5t 
and t' = OAt. For this value of t' the ground state 
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FIG. 1: (Color online) Temperature dependences of the in- 
verse antiferromagnetic susceptibility at t' /t = O.lt, U = 2.5t, 
and different values of the chemical potential (the list of the 
chemical potentials and fillings corresponds to the curves from 
top to bottom, the smallest ^ corresponds to upper curve). 
Dashed lines show the extrapolation of the inverse suscepti- 
bilities to the temperature region T < T^q" by polynomials 
of 6-th order 



FIG. 2: (Color online) Temperature dependences of the in- 
verse magnetic susceptibility at t' /t = O.lt, U — 2.5t and 
the incommensurate wavevector determined by a maximum 
T^, dashed lines show the extrapolation to T < Trg". Dot- 
dashed line shows the inverse commensurate susceptibility at 
^ = fi^^ ~ — 0.0375f. The inset shows temperature depen- 
dences of the inverse susceptibility with respect to d-wave 
superconducting pairing at different values of the chemical 
potential 



was previously found unstable with respect to antifer- 
romagnetic order and/or superconductivity at the fiU- 
ings close to van Hove band fiUing^t^ i^^i^^ . Tempera- 
ture dependences of the inverse antiferromagnetic sus- 
ceptibility (Q = (7r,7r)) obtained in the present approach 
for different chemical potentials are shown in Fig. 1. 
One can see that for large enough chemical potential 
fi > fi^^ « — 0.0375t = corresponds to van Hove 
band filling), the inverse antiferromagnetic susceptibil- 
ity monotonously decreases with decreasing temperature 
and vanishes at a certain temperature T^p. The value of 
TJp increases with increasing fi. 

Study of susceptibilities at the incommensurate wave 
vectors (see Fig. 2) shows that close to /i^^ (in the 
range -OMt < ^ < -0.02i) we have > T^p for 
some Q — (tt, tt — (5). Therefore an instability with re- 
spect to incommensurate, rather than a commensurate 
magnetic order is expected in this interval of fi. At 
/i = /i^ = —0.06t we obtain Tq = 0, which shows exis- 
tence of a quantum critical point below half filling. Near 
the quantum critical point we find Xc^ ~ T, which is sim- 
ilar to the result of the Hertz-Moriya-Millis theorjs^ii^. 

The behavior of the inverse susceptibility with respect 
to the d-wave superconducting order is shown in the in- 
set of Fig. 2. Similarly to the inverse antiferromagnetic 
susceptibility, it monotonously decreases upon lowering 
temperature, with a different temperature dependence. 



The obtained phase diagram is shown in Fig. 3 
and contains antiferromagnetic, incommensurate mag- 
netic and superconducting phases. Away from half filling 
the commensurate antiferromagnetic order is expected to 
be unstable towards phase separation (to hole-rich and 
hole pure regions) although this possibility can not 
be verified in the present approach. The obtained value 
of T^gQ monotonously increases with increasing density 
for n < 0.94. Deeper in the antiferromagnetic phase 
the superconducting transition temperature is somewhat 
suppressed. The origin of this suppression comes from 
the competition between antiferromagnetic and super- 
conducting fluctuations. The coexistence of supercon- 
ductivity and antiferromagnetism, which is possible in 
the interval 0.87 < n < 0.94, can not be verified in the 
present approach. 

The region of the incommensurate phase obtained in 
Fig. 3 is much narrower, than that expected in the mean- 
field approachesi^^, which predict incommensurate in- 
stability in the most part of the phase diagram. In fact, 
accurate mean-field investigationai^iiSiii^ show, that sub- 
stantial part of incommensurate state in the mean-field 
approach is unstable towards phase separation into com- 
mensurate and incommensurate regions and therefore 
qualitatively agree with the renormalization group ap- 
proach. The presence of incommensurate phases within 
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FIG. 3: (Color online) Phase diagram at t' /t — 0.1 and 
U = 2.5t. The temperatures of the crossover into regime with 
strong antiferromagnetic, incommensurate magnetic, and su- 
perconducting fluctuations are marked by squares, triangles 
and circles, respectively, PS denotes a possibility of phase 
separation of the antiferromagnetic case. Dashed line (stars) 
show the temperature Tao", at which the fRG flow is stoped 



the renormalization group approach was noticed previ- 
ously for t' = in Ref.^. 

The density dependence of Tafm ("■) and T^gc (") ' 
tained in Fig. 3, is similar to that of the antiferromag- 
netic and superconducting gap components in the elec- 
tronic spectrum, Aafm(?i) and Adsc(f^), recently ob- 
tained within the combination of functional renormal- 
ization group approach and mean-field theorylS. Slower 
decrease of T^gdn) when going into the antiferromag- 
netic phase in the present approach is explained by the 
fact that in the present approach magnetic and supercon- 
ducting fluctuations are weaker coupled in the absence 
of spontaneous symmetry breaking, since the latter leads 
to opening a gap in the electronic spectrum at the Fermi 
surface. Contrary to the study of RefJ^ we included in- 
commensurate phases in our analysis. 

At U = 3.5t we obtain similar behavior of the magnetic 
and superconducting susceptibilities near the quantum 
critical point; the resulting phase diagram is shown in 
Fig. 4. Compared to the case U = 2.5t, the phase di- 
agram has broader region of the incommensurate phase. 
The crossover temperature into regime with strong super- 
conducting fluctuations approximately follows that for 
the incommensurate fluctuations, implying that the su- 
perconductivity in this case is possibly caused by in- 
commensurate spin fluctuations. To clarify this point, 
we plot in Fig. 5 the momentum dependence of the 
superconducting gap, obtained from the Bethe-Salpeter 



FIG. 4: (Color online) Phase diagram at t' /t = 0.1 and U = 
3.5t. The notations are the same as in Fig. 3 

analysis^*'. We see that the shape of the gap, calculated 
for U — 3.5t shows stronger deviation from the d-wave 
form, than for U ~ 2.5t, which indicates possible role of 
the incommensurate fluctuations in this case. 



IV. CONCLUSION 



We have investigated temperature dependence of the 
commensurate and incommensurate magnetic suscepti- 
bilities, as well as the susceptibility with respect to the 
d-wave pairing in the fRG framework, which allowed us 
to obtain the phase diagrams of the Hubbard model at 
different U. We obtain an intermediate phase with strong 
incommensurate fluctuations between the commensurate 
and paramagnetic phases, the former is characterized by 
a wavevector Q = (tt, tt — J). The size of the incommensu- 
rate phase increases with increasing interaction strength. 
The tendency towards incommensurate order near mag- 
netic quantum phase transition comes from the absence 
of nesting of the Fermi surface at flnite t' . The corre- 
sponding profile of static noninteracting spin susceptibil- 
ity Xo(Q) is however almost flat near Q = (7r,7r) (see, 
e.g. Refj^i) showing that one can not restrict oneself 
to fluctuations with only one certain Q, as assumed in 
HMM theory. At the same time, the obtained size of 
the incommensurate phase is much narrower, than ob- 
tained in the mean-fleld approache o^^'^^ , which is ex- 
plained by existence of a phase separation in both ap- 
proaches. Near the quantum critical point the inverse 
magnetic susceptibility with respect to the preferable 
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FIG. 5: Angular dependence of the superconducting gap for 
U = 2M, n = 0.87 (dot-dashed hue) and U = 3.5t, n = 0.84 
(soUd hne), t' /t = 0.1. Dashed line shows the standard A = 
{coskx — cos ky)/ A dependence 

order parameter shows in fRG approach almost linear 
temperature dependence, similar to that in HMM the- 
ory. The electron-paramagnon interaction, not consid- 
ered in the present study, may however change the crit- 
ical behavior of the susceptibility. Note that recently 
the incommensurate magnetic fluctuations were also con- 
sidered within the renormalization-group approach for 
fermion-boson model^^, where similar results were ob- 
tained. 

While the Mermin- Wagner theorem states no sponta- 
neous breaking of continuos symmetry in two dimensions 



at finite T, we have obtained finite temperature of van- 
ishing inverse magnetic and superconducting susceptibil- 
ities, which is the consequence of the one-loop approx- 
imation, considered in Eqs. As we argue in the 
Introduction, the obtained temperatures T*^ should be 
considered as a crossover temperature to the regime with 
strong magnetic fluctuations and exponential increase of 
the correlation length. 

The patching scheme invoking the projection of the 
vertices to the Fermi surface, used in the present renor- 
malization group study, may have some influence on the 
phase diagram. We expect, however, that this influence 
does not modify the phase diagram strongly. This is con- 
firmed by the recent two-loop studyS^ which necessarily 
includes corrections to the effect of the projection of ver- 
tices and shows that the effects of these corrections and 
the two-loop corrections to large extent cancel each other. 

The non-analytical corrections to the susceptibility 
and electron-paramagnon interaction vertices may be- 
come important near quantum phase transitions^. These 
corrections are however expected to produce much 
weaker effect, than the effects of the band dispersion con- 
sidered in the present paper. Investigation of the role of 
these corrections in the presence of van Hove singularities 
has to be performed. 

Application of the method considered in the present 
paper to ferromagnetic instability and detail comparison 
of the results of the present approach with the mean-field 
approach and quasistatic approach of Ref.^'* also has to 
be performed. 
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